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Abstract 



The magnetoresistivity tensor of an interacting two-dimensional electron 
system with a lateral and unidirectional electric or magnetic modulation, in 
a perpendicular quantizing magnetic field, is calculated within the Kubo for- 
malism. The influence of the spin splitting of the Landau bands and of the 
density of states (DOS) on the internal structure of the Shubnikov-de Haas 
oscillations is analyzed. The Coulomb electron-electron interaction is respon- 
sible for strong screening and exchange effects and is taken into account in 
a screened Hartree-Fock approximation, in which the exchange contribution 
is calculated self-consistently with the DOS at the Fermi level. This approx- 
imation describes both the exchange enhancement of the spin splitting and 
the formation of compressible edge strips, unlike the simpler Hartree and 
Hartree-Fock approximations, which yield either the one or the other. 

71.45.Gm, 71.70.Di, 73.61.-r 



1 



Typeset using REVTpjX 



I. INTRODUCTION 



Modern techniques allow the fabrication of semiconductor heterostructures incorporating 
a two-dimensional electron gas (2DEG) and a lateral periodic electrostatic potential (electric 
modulation) and/or a periodic magnetic field (magnetic modulation). In the presence of an 
external, constant, and perpendicular magnetic field, the modulation lifts the degeneracy of 
the Landau levels. The resulting Landau-band structure determines various oscillations of 
the magnetoresistivities, which usually provide with the only accessible information on the 
modulation strength. 

There are at least three types of modulation effects on the magnetoresistivities. First, 
the Weiss commensurability oscillations in the quasi-classical regime of low magnetic fields, 
both for the electric,^ and, more recently for the magnetic modulations,!!! have attracted 
most of the experimental and theoretical work.iH Second, the peculiar subband structure 
generated by a two-dimensional superlattice, known as the Hofstadter butterfly, leads to 
another type of commensurability oscillations, inside the Shubnikov-de Haas (SdH) peaks, 
their observation being currently the aim of important efforts.! And third, at high magnetic 
fields, the profile of the density of states (DOS) associated to the energy dispersion of the 
Landau bands, together with the exchange-enhanced spin splitting, may also determine an 
internal structure of the SdH maxima. In the present paper we will discuss only this last 
type of effects. 

The experimental results have been obtained for the resistivity p xx of GaAs-AlGaAs in- 
terfaces with an electric modulation in the x direction, i.e. of a one-dimensional character. 
The first investigation was performed on a modulation created by holographic illumination.@ 
In the absence of the modulation the spin splitting of the second Landau level, with n — 1, 
could not be resolved in the magnetoresistivity, but with increasing the modulation ampli- 
tude, the evolution of a double-peak structure was observed, and it was attributed to the 
van Hove singularities (vHS) of the DOS. 

Further experiments, on higher mobility samples modulated by etching techniques, have 
clearly shown the spin splitting. For a weak modulation again a double-peak, but also a more 
complicated triple-peak structure have been detected in the spin-polarized SdH maxima of 
p xx corresponding to the band n = l.i For a stronger modulation the spin splitting vanishes, 
but the magnetoresistivity minima at low even filling factors change into maxima at higher 
even filling factors, simultaneously with the shift of the minima towards the odd filling 
factors.i This behavior has convincingly been explained by the cumulated effects of the 
overlapping vHS from adjacent Landau bands. In the experiments mentioned so far the 
modulation period has been in the range of 300-500 nm, and much larger than the magnetic 
length. 

Another series of recent measurements has been performed on modulated systems pro- 
duced by growing on vicinal surfaces.^ This technique generates an electric modulation of a 
much shorter period, about 30 nm. The resistivity p xx also displays a multi-peak structure 
which may be related to the vHS. The anisotropy of the resistivity, as well as an abrupt 
onset of the spin splitting for a high magnetic field, have been clearly evidenced.0 

In all these experiments a detailed interpretation of the results is difficult, and still 
insufficiently clear. Beyond the technological constrictions, the difficulties arise from the 
complicated relationship between the DOS and the magnetoresistivities, and from the non- 
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negligible electron-electron interaction effects. Therefore, a transport calculation combining 
on the same footing the electron-modulation, electron-electron and the electron-impurity 
interactions, and also thermal effects, is needed. 

The anisotropy of the modulation may result in a high anisotropy of the conductivity 
tensor. The conductivity a xx has a scattering component (inter Landau level), which depends 
quadratically on the DOS and vanishes in the absence of the impurities. The conductivity 
o yy has an additional band component (intra Landau level), determined by the dispersion of 
the one-particle energies, which diverges for weak electron-impurity collision broadening.ll! 
Therefore, in a high mobility system, the band conductivity may cover the DOS effects in 
a yy and in the related resistivity p xx X^ 

The Coulomb interaction also yields opposite effects. The tendency of the electrostatic 
screening is to reduce the energy dispersion imposed by the modulation on the effective 
single-particle states, and hence to increase the DOS. At the same time, the exchange inter- 
action lowers the energy of the occupied states, enhances the energy gaps, but also broadens 
the Landau bands, decreasing the DOS. We have recently calculated the energy spectr and 
the resistivity tensori for a modulated system, in the Hartree-Fock approximation (HFA) , 
and we could explain the spin splitting observed in the magnetotransport experiments. How- 
ever, other exchange effects have been overestimated in the standard HFA. 

The main artifact in the HFA results has been the appearance of strong short-range 
charge- density oscillations in the presence of a weak external modulation of period much 
longer than the magnetic length, for any filling factor. The reason is the competition of the 
Hartree interaction, of a repulsive character, with the Fock interaction, of an attractive char- 
acter, which in the presence of the modulation may excite high charge- density harmonics.^ 
This is reminiscent of the fact that the HFA predicts an instability of the homogeneous 2DEG 
against the formation of a charge density wave for any filling factor,0 whereas experiments 
indicate an inhomogeneous ground state (Wigner crystal) only for very low filling factors. 
Another consequence of the strong exchange energy is a strong exchange broadening of the 
Landau levels. A related implication is a substantial, qualitative contradiction between 
the HFA and the results of the HartreeBll or Thomas- Fermi calculations^ of edge states. 
While the latter predict compressible edge strips much wider than the magnetic length, but 
only the bare spin splitting (which is negligible for GaAsV the HFA gives considerably nar- 
rower compressible edge channels, but a strong splitting.EsFilli However, the experimental 
confirmation of wide edge channelsi! suggests the domination of the electrostatic effects. 

In order to avoid, or at least to minimize, these artificial features of the HFA, our 
previous attempt to include the Coulomb interaction in a magnetotransport calculation for 
a modulated systems has been limited to short modulation periods.! The steepness of the 
energy dispersion, on the magnetic length scale, can reduce the relative importance of the 
exchange interaction, such that the HFA may become qualitatively reasonable. 

In the present paper we want to extend our calculations to the situation when the mod- 
ulation period is much longer than the magnetic length. Therefore our efforts will be mainly 
focussed on the electron-electron interaction. Some preliminary results have already been 
reportedB Our approach is based on a screened HFA (SHFA), in which we include the influ- 
ence of screening on the exchange interaction. Although we consider only static screening, 
this leads already to the desired reduction of the exchange effects and avoids the artifacts 
of the bare HFA. The screening is mainly determined by the DOS at the Fermi levelB 
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Therefore, when the latter will be in an energy gap, the screening will be weak, such that 
the exchange enhancement of the Zeeman splitting will remain essentially like in the HFA. 
However, when the Fermi level will intersect a Landau band, the exchange effects will sub- 
stantially diminish. 

The Coulomb interaction is included in a transport calculation, within the standard Kubo 
formalism. For this purpose we need to consider an electron-impurity scattering mechanism. 
We will describe it within a self-consistent Born approximation (SCBA). 

The realization of a magnetic modulation with a period of a few hundred nanometers, and 
a large amplitude, of the order of 1 T is technically feasible.@ In the presence of a constant 
external magnetic field such a modulation will produce Landau bands and a charge-density 
response. Hence, screening and exchange effects will occur, as for an electric modulation. 
Even in the absence of relevant magetotransport measurements, we will include in our cal- 
culation such a periodic magnetic field. For simplicity we will discuss only the case when 
the modulation, electric or magnetic, is unidirectional and sinusoidal. 

The paper is organized as follows. In Sec. II we derive the self-consistent equations of the 
SHFA which give us the ground state of the system. Then, in Sec. Ill, we discuss the impurity 
scattering and the conductivity tensor. The numerical results of the transport calculation 
are presented in Sec. IV, and the conclusions are collected in Sec.V. Some technical details 
are given in two Appendices. 



II. SCREENED HARTREE-FOCK APPROXIMATION 

We combine the influence of the electron-electron and of the electron-impurity interac- 
tions on the single-particle states of the modulated 2DEG with the help of the average Green 
function, having the operatorial definition 

G(E) = (G-(E)) imp = E _ {H0 + l ee - Sei(E)) , (2.1) 

with the following notations: {■■■) imp stands for the average over all the impurity configura- 
tions; G^{E) = (E — H ± z0 + ) _1 , with H a generic one-body Hamiltonian of the interacting 
2DEG with impurities; H° is the Hamiltonian of the noninteracting 2DEG without impu- 
rities; S ee and S e * are the self-energy operators determined by the electron-electron and 
electron-impurity interactions respectively. 

In our case the noninteracting Hamiltonian has the form 

H° = -^-(p + eA(r)) 2 + V cosKx - %fx B B{x). (2.2) 
2m* 2 

The electrons are located in the plane {r = (x,y)}. B(x) is the projection of the magnetic 
field along the z axis and it may have a periodic component, B(x) = B + Bi cos Kx, similar 
to the periodic electrostatic potential, V(x) = V cosKx. We choose the vector potential in 
the Landau gauge, as imposed by the symmetry of our system, 

A(r) = (0 1 B x + -±sinKx). (2.3) 
K 
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We have also included in H° the Zeeman term, where a = + for spin-up and a = — for 
spin-down states, g is the bare, band structure g-factor, and \ib the Bohr magneton. 

We use the eigenf unctions of H° corresponding to the unmodulated system, i.e. for 
B\ = and V = 0, as the basis for the one-particle Hilbert space. These functions are the 
well-known Landau wave functions, f n x (x, y) = L~ l ' 2 exp{— iX y /l 2 )f nXo (x), where f n x ( x ) 
are the one- dimensional harmonic-oscillator wave functions, with n = 0,1,..., centered at 
the position Xq, called the center coordinate. Here I is the magnetic length, and we will 
denote by uj c the cyclotron frequency, both corresponding to the uniform component of the 
magnetic field B , I = (h/eBo) 1 ^ 2 and cu c = eB /m*. In order to simplify the notations 
we keep the same symbol for the two- variable wave function which depends on both spatial 
coordinates x and y, e. g. f nXo (x,y), and for the reduced wave function depending only 
on x, fnX Q {x). The distinction will be made by the number of variables specified inside the 
brackets. The plane- wave factor has been normalized to the macroscopic length L y . The 
matrix elements of H° are given in the Appendix A. 

We will assume randomly distributed impurities, such that the modulated system is 
invariant to translations along the y axis. Consequently, the dependence on y of the in- 
teracting, effective, one-particle wave functions, also factorizes in a simple plane wave. For 
these wave functions we use the notation ip n ax (x,y) = L y 1 l 2 exp(— iX y /l 2 )^ na - Xo (x), and 
we expand them in the Landau basis, 

tpnaX {x) =J2 C nn'( cr , X o)fn'X (x) . (2.4) 
n> 

For the interacting, unmodulated system, we have c nn >(<r, X ) = 5 nn >, and the dependence 
on the spin label arises as long as the exchange interaction and an external modulation are 
simultaneously present. The periodic fields broaden the degenerate Landau levels E na into 
energy bands, E no . XQ , which we find by solving the eigenvalue problem 

(H° + S ee ) Vw {x, y) = E nuXo ip nuXa (x, y) . (2.5) 

The averaged effect of the impurities consists in the spreading of the effective single- 
particle energies around the energy spectrum given by Eq.( |2.5| ), and thus in an additional 
energy broadening. The statistical weight of an arbitrary energy E is given by the spectral 
function, 

PnaXo(E) = -Im (ij naXo | G(E) I ij naXo ) , (2.6) 

7T 

and the contribution of the effective state (naXo) to the filling factor v can be defined as 

vnax,= J dE PnaXo {E)T{E) , (2.7) 

J-'(E) = [exp(E — n)/T + l] -1 being the Fermi function, with /i the chemical potential and 
T the temperature, such that 

^ f a dX 

f = £ / — *Wo- (2-8) 



5 



We will discuss the electron-impurity interaction model in the next section. Clearly, in the 
absence of the impurities p na x Q (E) = 5(E — E naXo )- 

For calculating the electron-electron self-energy we start with the form given by the 
HFA.@il In the Landau basis the matrix elements n r a , x i do not mix the spin and the 

center-coordinate quantum numbers, and can be written as 

Kn> 0, *o) = "rnrYo J dr dx' f* X() (x, y) V4rY {J, y')u(r - r) 

ttitYq 

X [fn'X ( X > y)4>mTY , u') x t) (x',y')i) mT Y (x,y)} 

= T%F(X ) + T%f(v,X ), (2.9) 

where u(r) = e 2 /(n | r |) is the Coulomb potential with k the dielectric constant of the 
semiconductor background. The two terms in the square brackets of Eq. fl2.9Q define the 
Hartree (direct) and the Fock (exchange) contributions, indicated by the superscripts H 
and F. The Hartree term can be rewritten in terms of the particle density, 

n ( X ) = U rnrY | IpniTYo(x) ?= J2 H P COS P Kx ' ( 2 ' 10 ) 

mrYo p>0 

the last form being the Fourier expansion appropriate to our external fields. The Fourier 
coefficients of the Hartree self-energy are given in Appendix A. 

Using Eq. fl2.4T ) and the Fourier transform of the Coulomb potential, w(q) = 2ne 2 / (k \ q |), 
the Fock self-energy becomes 

^nn^( CT > -^o) = — X! V moYo X! C m mi (& , ^o) c mm 2 ( a j ^o) 
mYo mrm 2 

X / ^f^fnx | e^ r I f miYo )(fm 2 Y \ e** \ f n >x ) ■ (2.11) 

The Fourier expansion of S^f(cr, X ) is also given in Appendix A. 

In order to overcome the artifacts of the HFA mentioned in the Introduction, we screen 
the Coulomb potential in the exchange self-energy, substituting w(q) in Eq.( [2.11|) by 

«(q) = --^r, (2.12) 
K qe(q) 

where e(q) is the static dielectric function. We thus neglect the dynamic screening effects. In 
the spirit of the random-phase approximation, e(g) = 1 — x(q) being the dielectric 

susceptibility given by the well-known Lindhard formula,^ to be evaluated self-consistently 
with the effective states resulting from Eq.( [2.5|) . In the diagrammatic picture, Fig.l, that 
means the series of the polarization loops — previously contained only in the Hartree part 
of E ee — is now included, via Eq. ([2.12|) , also in the Fock part. 

For a two-dimensional system in a perpendicular magnetic field one can identify two 
components of the static dielectric susceptibility, 

X(q) = Xi(q) + X2(q) , (2.13) 

corresponding to the intra- and to the inter- Landau level transitions respectively!!. In other 
words, Xi{q) describes the electrostatic response due to the electron redistribution around 
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the Fermi level, under the action of an electric field of an arbitrary wave vector q, while 
Xi{q) gives the response due to the distortion of the effective wave functions. 

For Xi we use t ne expression derived for homogeneous systems by Labbe0, which we 
average over the Brillouin zone determined by the external electric or magnetic modulation: 



a dX dv na x 
i dfj, 



F„ 



1 2 



(2.14) 



where we have used Eq.flAl]). We expect that this approximation is appropriate for weakly 
modulated systems characterized by IV <C ahu c , or IB\ <C aB . Within the same procedure 
we evaluate X2 as0 



f ,m 

± nn' I „ 



f a dX ® VncrX ~ ^n'aXp (2 15) 

a Jo a E na x — E n / a x 



For a sufficiently high magnetic field B , that is for sufficiently large energy gaps E n+1(T x — 
E na x , and for a small q, x^(o) is typically negligible with respect to Xi(<?)> except eventually 
when there are very few states at the Fermi level such that Xi(qO is vanishingly small. 
Since for our weak modulations, in the limit q — > 0, Xi(l) becomes proportional to the 
thermodynamic DOS, i.e. dn/dfi, we can say the screening of the exchange interaction is 
dominated by the intra-Landau level transitions, which at low temperatures are determined 
by the DOS at the Fermi level. 

This screened Hartree-Fock approximation (SHFA) has been used in the papers by Ando, 
Ohkawa and Uemura for explaining the exchange enhancement of the Zeeman0 and of the 
valley splittingil in homogeneous Si MOS systems. Further improvements have incorporated 
the dynamic screening, within the plasmon-pole approximation, and the Coulomb-hole cor- 
relation effects, for the calculation of the photoluminescence energy in n-doped GaAs quan- 
tum wellsllS. Nevertheless, in order to keep the consistent treatment of groundstate and 
transport properties tractable, we will neglect such additional corrections in the study of 
the modulated systems. In particular, even without dynamic screening effects, our results 
for the self-energy for a homogeneous 2DEG in a GaAs-AlGaAs interface are qualitatively 
similar to those shown in Refs. 33 and 34. In Fig. 2 we compare the first five spin splitted 
Landau levels vs. magnetic field, i. e. the Landau fan, in the standard (bare) HFA with the 
results given by the SHFA. The material parameters are m* = 0.067mo, g = 0.4, k = 12.7, 
and the carrier concentration is chosen such that uB — 10 T, v being the filling factor. 
Both the Landau and the spin gaps are enhanced by the exchange interaction, but in the 
SHFA the Fock self-energy is strongly dependent on the DOS at the Fermi level, D(Ep). 
Thus the exchange interaction is screened for a high D(Ep), i.e. for noninteger filling fac- 
tors, but it may become very large for an integer filling. Since the exchange interaction 
is negative, the screening mechanism leads to the deep cusps in Fig.2(b), also present in 
the dynamic-screening calculations.il'il Like in the HFA, the largest energy gaps occur for 
integer filling factors. The thermal energy is much smaller than the energy gaps, and no 
disorder broadening is assumed. Therefore the gaps at integer filling factors in Fig. 2(b) are 
smaller than those in Fig. 2 (a) only because of the screening involved by the response of the 
wave functions, described by X2- 

We believe that the energy spectra, which we obtained within the SHFA for the mod- 
ulated system, are much improved with respect to those in the HFA. In Fig. 3(a) we show 
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the SHEA of the Landau bands generated by an electric modulation of a period much larger 
than the magnetic length. The Zeeman splitting is enhanced in the states around the Fermi 
energy, with replica in each of the upper and lower Landau bands, as in the HFA.0 The 
improvement consists in the recovering of the pinning effect, resulting from electrostatic 
screening,^ on the spin-splitted energy bands near the Fermi level. The pinning effect 
defines both compressible (dispersionless) and spin polarized strips. In the HFA, for the 
modulation parameters of Fig.3(a), strong short-range oscillations of the effective energies 
and of the charge density, with typical periods of 3 — 5/, would occur. Those oscillations 
would disappear only if the energy dispersion imposed by the external modulation would 
dominate the interaction effects, i.e. for a sufficiently short modulation period and/or for a 
sufficiently large modulation amplitude.^ In sufficiently strong modulaton or steep confine- 
ment potentials the exchange enhancement of spin splitting is suppressed. With decreasing 
slope of the potential, the enhancement recovers suddenly, and a large spin splitting of the 
Landau bands occurs nearly symmetrically with respect to the Fermi energy in such a man- 
ner, that both spin polarized bands penetrate the Fermi energy with large slopes. Thus, the 
compressible strips obtained in the Hartree approximation are destroyed in the bare HFA 
and replaced with spin polarized incompressible strips. This spin polarization of the edge 
states may occur spontaneously even if the bare g-factor vanishes, and has been discussed 
as a type of phase transitionr9c3 We believe that the occurrence of this spontaneous spin 
polarization is an artifact of the unscreened HFA, since our SHFA yields (as long as a > I) 
only a very smooth and gradual reduction of the spin splitting with increasing modula- 
tion strength, with compressible spin polarized strips (energies pinned to the Fermi energy) 
instead of incompressible ones. 

In Fig. 3(b) we transpose the results shown in Fig.3(a) for a magnetic modulation. In 
this case the Hartree (direct) screening is weaker, due to the weaker, purely quantum- 
mechanically, coupling of the charge density to the nonuniform magnetic field, and conse- 
quently the pinning effect is less pronounced in Fig.3(b). The peculiar Hartree response 
to a nonuniform magnetic field, which becomes negligible in the classical limit (cyclotron 
radius at the Fermi energy larger than the modulation period) will be discussed in detail 
elsewhere.il In the present paper we restrict ourselves to a strong uniform component, B , 
such that the electrostatic and the exchange effects due to the periodic component are similar 
to those due to the electric modulation. 

III. CONDUCTIVITIES 

In order to calculate conductivities we have to consider the electron-impurity interaction. 
For randomly distributed impurities, the self-energy S ei is diagonal with respect to Xq, 
and obviously also with respect to o. Since we are mostly interested in the effects due 
to the electron-electron interaction, we simplify the calculation of S ei by employing the 
phenomenological ansatz in which one ignores the dependence of its matrix elements on 
the Landau quantum number n and on the center coordinate X , that is ^nax riaX (E) ~ 
Yf£(E)5 nn i . In the SCBA, Fig.l, we need to solve the equationlll 

V Jo a 
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ra dX 1 

= r y/ i (31) 

^Jo a E-E ncrXo -Xf(E) ' v ' ; 

We also take a simple parametrization of the electron-impurity interaction energy, T = 
1\J B [Tesla] [meV]. Using Eq. (|3.1| ) the DOS can be written in the form 

D(E) = ^ 2 D Y / ImK t (E), (3.2) 

where D = m* /2nh 2 is the DOS per spin level for the homogeneous system without mag- 
netic field. 

The standard Kubo formula for the conductivity tensor can be written as 

roo glut J . 

(Tap(w) = e J dt e- £t Tr(F(H)[fp, e^ Ht j a e^ H % mp , (3.3) 

where a, (3 — x, y, r a = a, and e — ► + . The operator r is given by 

r = ~[if,r] = v° + v ex , (3.4) 

where v° is the usual form of the velocity operator, resulting from the commutator of the 
noninteracting Hamiltonian H° with the position, 

v° = ^(p + eA(r)) . (3.5) 

All the other terms of the full Hamiltonian H commute with r except the Fock self-energy, 
due to its nonlocal character. According to Eqs.( |2.9| ) and ( |2.10| ), the Coulomb interaction, 
gives an exchange contribution to the velocity operator, of the form 

y e* = V E ee,F r ]_ (3.6) 

The matrix elements of the velocity operators ( |3.5| ) and ( |3.6| ), in the Landau basis, can be 
found in Appendix B. For the current- density operator we adopt the definition 

' r, (3.7) 



L X Ly 

and we will call the term corresponding to Eq. fl3.6p "exchange current". After a straightfor- 
ward manipulation the Kubo formula can be put in the form 

<7ap(u) = — [Kp{u) - A a/9 (0)] , (3.8) 

with 

o2 r oo 



— e r°° 
\ a p(u) = -r-r \ dET{E) 

±J X ±Jy J — OO 

x Tr{(6(E - H)r a G + (E + hu)fp) irnp + (r?G-{E - Kw)rJ{E - H)) imp ) . (3.9) 
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Of course r depends on the impurity configuration, via the exchange contribution. 
However, since in the following we will restrict ourselves to a weak disorder, the de- 
pendence of Eq.( |3.6| ) on the disorder broadening will be negligible. Therefore we ne- 
glect the impurity effects on r in Eq. (|3.9| ). Then, making use of the formal relation 
5{E — H) = [G~(E) — G + (E)]/2iri, we find that in the static limit, uj — > 0, we have to 
average combinations like (G ± (E)r a G ± (E))i mp . The result can be written as 

(G ± (E)r a G ± (E)) imp 

e / . i 



L X Ly 



(G^E))^ (f a + % -[Y, ei {E),r a ]\ {G ± {E)) impi (3.10) 



where we have used a Ward identity to express the current vertex correction due to the 
impurities with the help of the corresponding self-energy operator.il Since we have assumed 
T, ei (E) to be a trivial (spin-dependent) c- number, the impurity vertex correction vanishes. 
But we see that the Coulombian vertex correction is automatically included in Eq. (|3.10|) via 
the exchange current. 

The Hall conductivity can now be written as^I 

o xy = -a yx = — dETiE) 



o a 

nn'a 



xlm G naXo (E) Re -^G n , aXo (E) . (3.11) 

Since we treat the disorder in the SCBA, we cannot consider localization effects, and hence 
we cannot expect real Hall plateaus. Moreover, if the disorder broadening of the Landau 
levels is much smaller than the energy bands, Eq.( |3.11| ) can be well approximated by the 



limit of a vanishing electron-impurity interaction, 7 — > 0. In this limit a xy takes the more 
familiar form 

ihe 2 f a dX _ 

(4>naX I X I 1pn'aX )(tpn'aX \ V \ 4>naX ) 
{-&naX — &n'aX ) 

in which the energy gaps are explicitly evidenced. 

It is instructive to mention that neglecting the exchange current one obtains deviations 
from the well known quantized values, a xy = {e 2 /h) x (integer), for an integer filling fac- 
tor. Let us consider the simplest case, with no external modulation. In the absence of the 
Coulomb interaction the energy gaps are determined by the cyclotron energy. The cyclotron 
frequency squared, in the denominator of Eq. (|3.12|) , is compensated by the cyclotron fre- 
quencies introduced in the numerator by the velocity matrix elements, Eq. (|B~T|) , and for a 



sufficiently low temperature one obtains the simple Drude formula, a xy = (e jh)v. When 
the Coulomb interaction is present, the energy gaps are enhanced, but a similar enhance- 
ment occurs in the numerator, determined by the exchange term of the velocities, Eqs.( [3.6| ) 
and (|B3|) , and the Drude formula remains valid. Obviously, when the exchange interaction 
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is screened the exchange currents may be small. In our calculations, for the modulated sys- 
tems, their contribution to the conductivities have been of the order 30-50 % in the HFA,1 
but of no more than 10 % in the present SHFA. 

The longitudinal conductivities can be transformed from Eqs. (|3.8|) - (|3.9|) to 

a aa = J dE (^-^jj VaaiE) , 

a aa (E) = ^— [ — J2 I (VW I r a I ipn'aXo) 1 2 
L^7l z Jo a i 

nn'a 

x Im G naXo (E) Im G n > aXo (E) , (3.13) 

The relation between the longitudinal conductivities and the DOS is complicated. Both 
the diagonal conductivities a xx and a yy have inter-level components, corresponding ton ^ n' 
in Eq.( [3.13j ), also known as scattering conductivities. One can show, by inspecting Eq.( [3.13j ), 



that for a weak disorder, 7 — > 0, the scattering conductivities become proportional to 
(TD(Ef)) 2 . Hence, they qualitatively reproduce the DOS profile, with van Hove singularities 
corresponding to the edges of the one-dimensional Landau bands. 

Due to the anisotropy of the system, an intra- level component of Eq.( |3.13D , with n = n', 
is nonzero for a yy only. It is called band conduct ivity,B0 being directly related to the 
dispersion of the Landau bands which yields equilibrium Hall currents. If the Hamiltonian 
is a local operator, the Hellman-Feynman theorem leads to 



dE ncrXo 
dX n 



-mu c (ij naXo I V I ^PnaXo) , (3.14) 



which is no longer valid when the exchange interaction is considered. However, if the latter 
is screened, the deviations from Eq.( |3.14j ) are not very important, and the band conductivity 
can still be understood in terms of the energy dispersion. Classically, it corresponds to the 
quasi-free drift of the centers of the cyclotron orbits along the y axis, parallel to the periodic 
electric field. Thus the band conductivity behaves contrary to the scattering components: 
It vanishes at the band edges, has maxima at the band centers, and diverges for a small 
disorder, like 7 -2 , becoming the dominant contribution to o yy . 

The structure of the longitudinal conductivities, a xx (E) = a^ ttermB (E) and a yy (E) = 
a scattenng ^ ^.band^^ £ Qr an i so i a ted, sinusoidal Landau band are qualitatively displayed 
in Fig. 4. In our calculations the two scattering conductivities are nearly equal. In Fig. 4(a) we 
assume a small band conductivity, such that both a xx and a yy show the two van Hove peaks. 
Reducing the disorder broadening the scattering components decrease, but at the same time 
the band conductivity increases, and a triple-peak structure evolves in a yy , Fig.4(b), while 
the shape of a xx remains unchanged. For an even smaller disorder a yy becomes dominated 
by the band conductivity, the central peak covering the vHS, like in Fig. 4(c). 



IV. DISCUSSION OF THE NUMERICAL RESULTS 

We perform the calculations of the effective electronic states within a numerical iterative 
scheme, starting from the noninteracting solution, and assuming a fixed number of particles 
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for determining the chemical potential. At each step we diagonalize the Hamiltonian matrix, 
Eqs.Q, (HI) and ([AToD , in L points of the half Brillouin zone, < Xq < tt/K, and then 
we compute the first L Fourier harmonics for Eqs.( [A4| ) and ( |A5| ) . We have taken L in the 
range 10-40, and we have mixed 5-10 Landau levels. Then, we use Eqs. ( |3.12|) and (|3.13|) to 
calculate the conductivity tensor. 

In order to resolve both the van Hove peaks in the modulated system and the spin 
splitting of the Landau bands, we need a very small disorder parameter. In Fig. 5 we show the 
magnetoconductivity tensor for a magnetic field B varied such that the Fermi level traverses 
the Landau bands with n = 1 and n = 2. In Fig. 5 (a) we consider a pure electric modulation 
of amplitude V — 15 meV and period a = 500 nm. For a better understanding of the 
conductivity oscillations, three qualitatively different energy band structures are indicated 
in Fig. 6, in a half Brillouin zone. The bands corresponding to opposite spin directions are 
separated for filling factors below 4, i.e. for B > 2.5 T, due to the exchange enhancement, 
and they partially overlap for higher filling factors. 

For the electric modulation the energy dispersion is strongly dependent on D(E F ), due 
to the strong screening contained both in Yf e,H and in S ee,F . All the Landau bands shrink 
whenever a band is intersected by the Fermi level. Therefore a full spin splitting — i.e. 
non-overlapping bands with the same n, but with different spin directions — can result even 
for a modulation amplitude much larger than Tioj c . 

Due to the small band width, and also due to the pinning of the band edges to the Fermi 
level, the resolution of the vHS peaks requires a very small disorder broadening. For our 
parameters the vHS can be distinguished in a xx , Fig.5, and their resolution even improves 
with decreasing magnetic field. The reason is the increase of the exchange interaction, self- 
consistently with the reduction of the DOS: D(Ep) decreases, the screening diminishes, 
and the exchange broadening is larger. Nevertheless, the screening may again increase for 
overlapping bands, like in Fig. 6(c). In that case two bands contribute to the pinning effect, 
and therefore the resolution of the vHS in the states (2, +, 0) and (2, — , tt/K), which yield the 
two maxima of a xx in Fig.5(a) for 1.9 T < Bo < 2.1 T, is poorer than in the states (2, +, tt/K) 
and (2, -, 0), that is for 2.4 T < B < 2.5 T and 1.6 T < B < 1.7 T, respectively. 

While the scattering component of a yy is nearly identical to a xx , the very small disorder 
parameter makes the band conductivity large. In Fig.5 ( St) (J yy IS thus in the situation depicted 
in Fig.4(c), with the vHS profile hidden by the band-conductivity single peak. Decreasing 
the magnetic field, due to the stronger energy dispersion the band conductivity further 
increases, contrary to the scattering term, which becomes negligible. For the overlapping 
bands of Fig. 6(c) the minimum of a yy at B = 2 T corresponds to the strongest screening, 
i.e. to the highest D(E F ). 

In Fig. 5(b) we show similar results obtained for a purely magnetic modulation which 
produces Landau bands equivalent to those of Fig. 6. We have chosen B\ = 0.25 T, the 
other parameters of the calculation being the same as those for Fig.5 (a). There are two, but 
unessential differences with respect to the electric modulation of Fig. 5(a). First, the energy 
dispersion increases with increasing Landau quantum number n. In a crude approximation, 
i.e. neglecting the Coulombian effects, as long as Kl -C 1 the energy bands are determined 
by the local cyclotron energy, E na x = ^{B$ + B\ cos KX )(n + |). Second, as mentioned 
at the end of Sec. II, the Hartree (direct) screening is weaker than in the electric case, and 
hence the strength of the exchange with respect to the Hartree interaction may increase 
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to some extent. Going through Fig. 5(b) from high to low values of the constant field Bq, 
the increase of the energy dispersion at the Fermi level due to the increase of the quantum 
number n is thus amplified by the exchange broadening. For filling factors higher than 4, 
four Landau bands overlap, i. e. those with n = 2 and n = 3. The order of the vHS 
observed in a XXl for B < 2.5 T is : (2, +, tt/K) (B = 2.48 T), (2, -, ir/K) (B = 2.10 T), 
(2, +, 0) (B = 1.84 T), (3, +, ir/K) (B = 1.72 T). The small scattering conductivity a xx for 
2.15 T < B < 2.40 T also reflects a strong energy dispersion. Additionally, the screening- 
exchange balance may generate weak DOS fluctuations near the band edges,! which can be 
seen as the shoulders of the band conductivity (via the Green function squared in the Kubo 
formula), in a yy for 2.0 T < B < 2.5 T. For lower magnetic fields, when D(Ep) increases 
due to the overlapping bands, the exchange effects are again small, and the maxima of the 
scattering conductivity corresponds to the minima of the band conductivity and vice-versa. 

For comparison with experiment we need to invert the conductivity tensor into resistiv- 
ities, 

_ ^VV _ ®xx _ @yx / , , \ 

Pxx — -p - , Pyy — -p - ) Pxy = i l 4 -^ 

where V = a xx a yy + a^y- m our calculations a xx a yy <C a^ y , and hence we have p X x, yy ~ 
a yy,xx/^ x y In Fig-5(c) we see that the resistivity measured perpendicular to the modulation 
reflects in fact the band conductivity, while the scattering conductivity alone can be observed 
only in the resistivity measured parallel to the modulation. 

In the experiments on the systems with an etched electric modulation, of parameters 
comparable to our model, more complicated structures of the resistivity p xx have been 
observedil The SdH maxima corresponding to a certain spin direction may show two and 
even three internal peaks. One possible explanation could consist in the more complicated 
exchange effects on the band conductivity, like the shoulders we have found for the magnetic 
modulation, Fig. 5(b). We have also obtained such effects for an electric modulation, with a 
shorter period, a = 100 nm, when the screening is considerably reduced]! Another possible 
explanation, which is supported by the present results of the SHFA, is that in the real system 
the band and the scattering conductivities may be comparable in magnitude, such that their 
superposition in a yy may lead to one to three peaks per Landau level, as illustrated in Fig.4. 
However, within our approximations, the price for the resolution of both the spin splitting 
and the vHS is a relatively large band conductivity. 

To our knowledge, a simultaneous measurement of both p xx and p yy , capable of indicat- 
ing the real magnitude of the two components of the longitudinal conductivities, has been 
reported, at high magnetic fields, only for the short-period modulated systems^ and p xx 
has been found considerably larger than p yy , but without a clear DOS structure. For short 
periods, the Hartree screening gets weaker and the Landau bands become wider, such that 
in principle we could cover several situations, including those of Fig.4, by varying the modu- 
lation amplitude and the disorder parameter 7. We cannot extend our SHFA to that regime, 
because of the steepness of the energy dispersion which is not compatible with the assump- 
tion about the quasi-homogeneous screening of the exchange interaction. Nevertheless, the 
results of the HFA may be satisfactory! 

We want now to discuss the situation of an electric modulation with a larger amplitude, 
such that the Landau bands overlap. In the experiment by Weiss et. al.i the double-peak 
structure of the resistivity p xx has been found, while the absence of the spin splitting for the 
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unmodulated system may be attributed to the combined disorder and thermal effects. We 
have recently confirmed such a possibility, within the SHEA, by choosing such a disorder 
parameter and temperature that both the exchange enhancement of the Zeeman splitting 
and the band conductivity have been suppressed, and p xx have been obtained as resulting 
from Fig. 4 (a) .0 

In the example of Fig. 7 we consider the situation in which neither the spin splitting nor 
the vHS are resolved in the low-energy Landau bands, but the vHS develop in the higher, 
overlapped bands. The SdH minima at filling factors around 2 and 4 are slightly shifted to 
higher magnetic fields. The energy gaps are small, nearly vanishing, but the adjacent bands 
are separated due to the pinning effect. For weaker magnetic fields, at even filling factors 
the vHS from the top of the Landau band below the Fermi level overlaps with that from 
the bottom of the upper band, resulting in the maxima of a xx for v — 8 and v = 10 and 
in the minima for the odd filling factors v = 7 and v = 9. The transition occurs around 
v = 6. Since the DOS decreases when the magnetic field is lowered, the screening becomes 
less efficient and the band width may increase when the Fermi level is in a band center, 
leading to a large band conductivity, wherefrom the maxima of a yy for v = 7 and v = 9. 
The switching from even to odd filling factors at the resistivity minima, which we show 
here only for p yy (cr xx ), has been recently observed for p xx and supported by a transport 
calculation similar to the present one, but with the Coulomb interaction neglected.! 

Finally we want to mention that the plateaus we obtain for the Hall resistivity, Fig. 7(b), 
do not correspond to plateaus of the Hall conductivity, since localization effects do not exist 
in our model, but to the wide minima of the longitudinal conductivities. In the noninter- 
acting approximation, and also in the HFA, for a fixed number of particles, the Fermi level 
jumps abruptly from one Landau band to the other. In the SHFA the exchange interaction 
is very sensitive to the variation of D(Ep), and since they are self-consistently determined, 
the stability of the chemical potential in an energy gap considerably improves. 

V. CONCLUSIONS AND FINAL REMARKS 

We have calculated energy spectra of a two-dimensional electron gas in an electric or 
a magnetic superlattice, and in a perpendicular magnetic field, beyond the standard HFA. 
The superlattice potential is smooth and the uniform magnetic field is strong, such that 
the electrostatic screening is very important. The essential element of our approach is the 
inclusion of screening in the exchange term, self-consistently with the DOS at the Fermi level. 
The numerical calculations are complicated and time consuming, such that approximations 
in treating the screening are unavoidable: we have only considered the static screening, in 
a manner appropriate to a quasi-homogeneous system. Improvements of our procedure are 
possible, e.g. along the lines of Refs. 32, 33, but even within the present approach the 
results drastically change with respect to the HFA. The obtained screening of the exchange 
interaction may be even somewhat too strong, since we found that our SHFA yields no 
charge- density- wave instability of the homogeneous 2DEG, even at very low filling factors. 

Using a large-amplitude modulation as model for edge states, we have shown that our 
SHFA results interpolate between the contradictory Hartree and Hartree-Fock approxima- 
tions: we obtained both compressible edge strips much wider than the magnetic length, and 
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an enhanced spin splitting. To the best of our knowledge this expected result has never been 
obtained within a microscopic many-body theory before. 

For a weak modulation, the standard HFA yields strong short-range oscillations of the 
Landau bands and the particle density, which have never been observed experimentally. Such 
oscillations completely disappear in our SHFA, and we believe that these present results are 
much more realistic. Clearly, we cannot expect a quantitative agreement with experiments 
since, besides the approximations regarding the electron-electron interaction, we have used 
a rather crude simplification of the electron-impurity scattering in our transport calculation, 
such that the disorder vertex correction vanishes. But qualitatively, we have obtained smooth 
internal structures of the SdH peaks determined by the interplay of the scattering and 
band conductivities, the ingredients of which have been carefully analyzed. These smooth 
internal structures compare much more favorably with the experimental results than the 
sharp structures near the Landau band edges obtained in the bare HFA.! 
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APPENDIX A: HAMILTONIAN MATRIX ELEMENTS 



With the notations z = (Kl) 2 /2 and 



/■;„,(:) = f ^ ) zW 2 e-*'*L»r n '(z), (Al) 



where L™, n (z) are the Laguerre polynomials,^ n, n' = 0, 1, and F nn i = (—1)" n ' F n > n , 
can write the matrix elements of the noninteracting Hamiltonian in the form 



we 



H° n ,(a, Xq) = 5 nn ,hu c 



1 1 (B x 
n H 1 — 

2 8z \B 



2 1 



a 



Sun 1 —9^bBq 



+ 



- ^F nn >(z) + \frm'F n -x^^{z) - \J(n + \){n' + i)F n+ltn >+i(z)j 



a 



V[V--gii B B 1 ) F nnl {z) 



cos (kXq + (n- nO^) 



hu c ( Bi 



8z 



F nn ,(Az) cos (2KX + (n-n')^ 



(A2) 



Due to the reflection symmetry of the external fields, we have the following parity rule, 



H nn ' (cr, X 



'-l) n n H° n ,(a, -Xq) 



(A3) 



which also holds for the Coulomb self-energy matrix elements, Eq. ( |2.9|) , which depend on 
the self-consistent wave functions given by Eq.( p.4|) . In order to write those matrix elements 
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in a form suitable for a numerical calculation, we expand the mixing coefficients in Fourier 
series, 

C n „/((T, X ) = E lrm>{v,P) COS (pKX + (fl - Tl')^ J , (A4) 

where we have taken into account the parity rule. 
The effective energies can thus be expanded as 

EnaXo = E e nap COS pKX , (A5) 
p>0 

and similarly the occupation numbers 



VnaX = E f™P COS P KX ■ (A6) 

p>0 



We define 

4 



j\n 1 n2n- i n A 
P1P2P3P4 



-/ dx Y[cos(piX + rii-). (A7) 

7T J-n i=1 Z 



With this notation, and using Eqs. flA4|) and ( |A6| ) , the Fourier coefficients of the particle 
density can be written, by inverting Eq. (|2.10|) , as 



2 — 5p0 y—y 

1U/ "' m 1 >0,cr = ± 
P1>0 

x E ^rT 3 ' '"r i " m3F - 2 ™3(P 2 ^)7 mim2 (cr,p 2 )7 mim 3(^P3). (A8) 

m2,m3>0 
p2,p3>0 

The matrix elements of the Hartree term, defined by Eqs.( |2.9| ) and ( |2.1U| ), can now be 
put in the form 

Kf(*o) = E ^F nn ,{p 2 z) cos (pKX + (n- n')|) , (A9) 

where we have used the Fourier transform of the Coulomb potential. The corresponding 
Fourier series for the Fock term of the Hamiltonian, Eq. (|2.11|) , can be found by directly 
searching for the Fourier amplitudes, and after a lengthy, but straight forward calculation, 
one gets 

x o) = ~ttK] E( 2 ~ W cos (pKXo + {n- n')^-) £ /m l0Pl 

P1>0 

X 



E ^rTrXr m2 P r" m3 ^^-3(p^)7 mim2 (a,p 2 ) 7mim 3(^P3) , (mo) 

m2,m3>0 
p2,p3>0 

where S denote the exchange integrals 
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. { mi\m 2 \\ 1/2 f°° , _,, 

S mini ,m 2 n 2 (t) = : — r / dqu{qV2) 

\ ni\n 2 \ ) Jo 

xe- g \ n i- mi+n2 - m > +1 J ni ~ mi -n 2+ rn 2 (tqV2^ (All) 

with J n (x) the Bessel functions. For the homogeneous system the self-consistent wave func- 
tions are identical with the Landau wave functions, so that r y ni n 2 ( 'yP) — ^nm^^poi an d the 
exchange integrals reduce to S , minii?nini (O).0@ 



APPENDIX B: VELOCITY MATRIX ELEMENTS 

The matrix elements of the noninteracting components of the velocity operators, Eq.(|3T 
in the absence of a magnetic modulation, are 

(v° a ) nn ,(X ) = vJ-^= (Vn+ l<5 n >+i + rfcVn' + 15„,„'+i) (Bl) 

with the notation (rj x ,r) y ) = (i, 1). When the magnetic field has a periodic component in 
the x direction, one has to add to v® the extra term 

Jj ' Bl -F nn ,{z) sin (kX +(n- n')^) . (B2) 



The matrix elements of the exchange current can be transformed in Fourier series 
like those of the exchange-interaction operator, as sketched in Appendix A. According 
to Eq.fl3.6p, the form of v^j,(cr, Xq) is similar to that of £^f(<r, X ), with the replace- 
ment u(r — r') — > (r — r')w(r — r') in Eq. (|2.9| ), and, after an integration by parts, with 
w(q) — > V q tt(q) in Eq. (|2.11| ), one obtains the following Fourier expansion: 

C( ff i x o) = 7= - E( 2 - <W sin (pKX + (n - n')£) E /Wow 
Hy/Znn p > V ^/ mi > 



x 



E ^T^r""™ 3 ^^ , (B3) 



m2,m3>0 
p2,p3>0 



in which 

1/2 



mini ,m2U2 



ni!n 2 ! / 



du(qy/2) 



dq 

xe 9 g^ 1 mi+ ™ 2 m-2+1 [J ni _ mi _ n2+m2 _ 1 (tg^/2) — J„, 1 _ mi _ ri , 2 + m2 +i(tq , V / 2)] 

x^r mi (? 2 )^r W2 (? 2 )- ( B4 ) 
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FIGURES 



FIG. 1. Diagrammatic representation of the combined screened Hartree-Fock and self-consistent 
Born approximations. The wavy lines are the Coulomb interaction and the dashed line is the 
electron-impurity interaction. 

FIG. 2. The Landau fan (a) in the HFA and (b) in the SHFA for a homogeneous GaAs system. 
The dashed line shows the chemical potential. The temperature T = 1 K. 

FIG. 3. Compressible and incompressible strips in the SHFA, for: (a) an electric modulation 
having a = 1000 nm and V = 200 meV with Bq = 6 T, and (b) a magnetic modulation having 
a = 500 nm and B\ = 1.2 T and Bq = 4 T. The Fermi level is indicated by the horizontal dashed 
line. T = 1 K. 

FIG. 4. Three possible profiles of the magnetoconductivity corresponding to a sinusoidal-like 
isolated Landau band. The full line shows a yy , the dashed lines a xx and the scattering component 
of cr yy , and the dotted lines the band conductivity. 

FIG. 5. The conductivity tensor calculated for (a) an electric modulation with a = 500 nm and 
V = 15 meV, (b) a magnetic modulation with a = 500 nm and B\ = 0.25 T, and (c) the resistivities 
corresponding to (a). The xx components are magnified by a factor of two. The temperature T = 1 
K and the disorder parameter 7 = 0.06. 

FIG. 6. Three energy spectra corresponding to Fig.5a. 

FIG. 7. (a) The conductivities, and (b) the resistivities, for an electric modulation with a = 500 
nm and V = 20 meV. The Hall components are reduced by a factor of two. The numbers inside 
the plot indicate the filling factors. T = 1 K and 7 = 0.20. 
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